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Abstract 

We study a system of self-gravitating massive fermions in the framework of the 
general-relativistic Thomas-Fermi model. We postulate the free energy functional and 
show that its extremization is equivalent to solving the Einstein’s field equations. A self- 
gravitating fermion gas we then describe by a set of Thomas-Fermi type self-consistency 
equations. 

Thermodynamical properties of the self-gravitating fermion gas have been extensively 
stndied in the framework of the Thomas-Fermi model [1-6]. The system was investigated in 
the nonrelativistic Newtonian limit. The canonical and grand-canonical ensembles for snch 
a system have been shown to have a nontrivial thermodynamical limit |^, . Under certain 

conditions this system will nndergo a phase transition that is accompanied by a gravitational 
collapse 1^, ^ which may have important astrophysical and cosmological implications [^ . 

In this paper we formnlate the general-relativistic version of the model. The effects 
of general relativity become important if the total rest-mass of the system is close to the 
Oppenheimer-Volkoff limit [0. There are three main featnres that distingnish the relativis¬ 
tic Thomas-Fermi theory from the Newtonian one: i) the eqnation of state is relativistic 
ii) the temperatnre and chemical potential are metric dependent local qnantities Hi) the 
gravitational potential satishes Einstein’s held eqnations (instead of Poisson’s eqnation). 

Let ns hrst discnss the general properties of a canonical, self-gravitating relativistic hnid. 
Consider a nonrotating hnid consisting of N particles in a spherical volnme of radius R in 


1 


equilibrium at non-zero temperature. We denote by , p, p, n and a the velocity, pressure, 
energy density, particle number density and entropy density of the fluid. Following Gibbons 
and Hawking [Q we postulate quite generally the free energy of the canonical ensemble as 

F = M - (1) 

where M is the total mass as measured from infinity, S is a spacelike hypersurface that 
contains the fluid and is orthogonal to the time-translation Killing vector held which is 
related to the velocity of the huid 

F = (2) 

The entropy density of a relativistic huid may be expressed as 

a=^{p + p-pn), (3) 

and T and p are the local temperature and chemical potential. The metric generated by the 
mass distribution is static, spherically symmetric and asymptotically hat, i.e. 

ds^ = — X^dr‘^ — r‘^{d9‘^ + sinOdcj)'^). (4) 


^ and A may be represented in terms of the gravitational potential and mass 

e 


A= 1- 


2m(r) 


- 1/2 


with 


nr 

m{r) = / dr' A'Kr'‘^p{r') . 
Jo 


The equation of hydrostatic equilibrium P, 10 


duP=-{p + p)^ 

together with the thermodynamic identity 

dt = Tii + tdl, 

n 1 n 1 


(5) 

( 6 ) 

(7) 

( 8 ) 

( 9 ) 
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implies 


rj = To ; = /ici, 


( 10 ) 


where Tq and po are constants eqnal to the temperature and chemical potential at inhnity. 
To may be chosen arbitrarily, but /xq is an implicit functional of ^ due to the constraint that 
the total number of particles 

( 11 ) 


= N 


is hxed. 

Based on Eq. (p!0D the free energy may be written in the form analogous to ordinary 
thermodynamics 

F = M -TqS (12) 

with M = m{R) and the total entropy S dehned as 


S = f drAwr^X^ip + p) - 

JO 1 In 


(13) 


where we have employed the spherical symmetry to replace the proper volume integral as 


cR 




drAnr^X. 


(14) 


The following theorem shows that the quantity dehned by Eq. (0) is indeed the free 
energy of a self-gravitating relativistic huid at hnite temperature. 


Theorem 1 Among all momentarily static, spherically symmetric configurations {^(r), m(r)} 
which for given temperature Tq at infinity contain a specified number of particles 



Anr'^dr A(r)n(r) 


N 


(15) 


within a spherical volume of given radius R, those and only those configurations that extrem- 
ize the quantity F defined by (0) will satisfy the Einstein’s field equation 


dr 




m + Airr^p 
r{r — 2m) ’ 


( 16 ) 


with the boundary condition 


2M 


1/2 





Proof. From Eqs. (|^ and (|^ by making nse of the identity (|^), and the fact that 
= 5{^q/Tq) and that N is hxed by the constraint (|^, we hnd 


6F = 6M — [ dr 4:7ir‘^^{p + p)S\ — [ dr 4:'nr‘^\-;^5p. (18) 

Jo T Jo T 

The variations 6X and 6p can be expressed in terms of the variation 6m{r) and its derivative 

d6m 


This gives 


6F = 6M- 


rR 


dr 

.Tn 


= Anr^Sp. 


dr47Tr^:!^{p + p)^6m- f^drx'^^. 

Jo T dr 


By partial integration of the last term and replacing Tq/T hj f we find 


6F=[1- X(R)f{R)] 6M- , 

Jo 


dr 


. 9^/ ^ dX d 

p + p -- 

am dr 


6m, 


(19) 

( 20 ) 

( 21 ) 


where 6m{r) is an arbitrary variation on the interval [0, R], except for the constraint Sm{0) = 
0. Therefore 6F will vanish if and only if 


Airr^fip + p)^ - = 0 

am dr 


( 22 ) 


and 


1 - X{R)f{R) = 0. (23) 

Using Eqs. (|§) and (|^) we can write Eq. (p2| ) in the form (pTl), and Eq. ( p3[) gives the 


desired bonndary condition ([T7|). Thus, 6F = 0 if and only if a configuration {f,m} satisfies 
Eq. ([T^ ) with (pTD as was to be shown. 

Remark 1. A solutions to Eq. ([I6|) is dynamically stable if the free energy assumes a 
minimum. 

Remark 2. Our Theorem 1 is a hnite temperature generalization of the result obtained for 
cold, catalyzed matter [|^ . 

We now proceed to the formulation of the general-relativistic Thomas-Fermi model. Con¬ 
sider the case of a self-gravitating gas consisting of N fermions with the mass mf contained 
in a sphere of radius R. Given the temperature at infinity Tq, the following set of self- 
consistency equations defines the Thomas-Fermi equation: 

1 Equation of state 

/■°° dJq 1 , 

""^^Jo (^ 1 + 
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(25) 


P = 9 


d?q 


E 


p = gT 


0 (27r)^ 1 + ’ 

ln(l + , 


roo (]3q 


(2vr) 


(26) 


where g denotes the spin degeneracy factor, T and g, are local temperatnre and chemical 

)^), and E = 

m + dvrr^p 


potential, respectively, as dehned in Eq. ( 0 ). and E = ^mj + 

2 Field equations 


_ . 

dr r(r — 2m) ’ 

dm , n 
—— = 4:7ir p, 
dr 


with the bonndary conditions 


?(fl) = 1 


3 Particle number constraint 


rR 

/ dr 47rr^(l — 2 m/r)“^Z 2 

Jo 


2m(/2)y^^ 


R 


m{0) = 0. 


(27) 

(28) 

(29) 

(30) 


One additional important reqnirement is that a solntion of the self-consistency eqnations 
(^4[) to (^^ have to minimize the free energy dehned by Eq. (|^ 

Next we show that, in the Newtonian limit, we recover the nsnal Thomas-Fermi eqnation. 
Using the nonrelativistic chemical potential pnr = go — RJf and the approximation ,^ = ~ 

1 + ip, E mf + q‘^/2mf and m/r <C 1 , we hnd the nsnal Thomas-Fermi self-consistency 
equations |^, || 


P 

n = -= q 

mf 


'fl h + exp) ^ 


'0 (27r)3 


, mf Pnr, 

- -j- - ip — - 

2m/To To To 


-1 


dip m dm , , 

-r = -}- = dvrr p, 


dr 


dr 


t(^) = -^^^; m(0) = 0, 

eR 

I dr inr^nlr) = N. 

Jii 


(31) 

(32) 

(33) 

(34) 
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The free energy (|^ in the Newtonian limit yields 

1 rR 


F = 


1 

irifN + unrN -/ dr Anr'^nip — / drAnr'^p 

2 Jo Jo 


with 




(27r 


2m/To 


rrif 

—-{ 
To 


dr Anr'^p 

(35) 

+ 7 )) . 

(36) 


which, np to a constant, eqnals the Thomas-Fermi free energy [^. 

A straightforward thermodynamic limit iV —> cx) as discussed by Hertel, Thirring and 
Narnhofer p|, ^ is in our case not directly applicable. First, in contrast to the non-relativistic 
case, there exist, at zero temperature, a limiting configuration with maximal M and N 
(Oppenheimer-Volkoff limit) and something similar can be expected at T 7 ^ 0. Second, the 
scaling properties of the relativistic Thomas-Fermi equation are quite distinct from the non- 
relativistic one. The following scaling property can be easily shown: 

If the configuration {^(r),m(r)} is a solution to the self consistency equations (^) to 
(|30[), then the configuration = ^{A~^r),rh = Am{A~^r)-, A > 0} is also a solution 
with the rescaled fermion number N = A^^'^N, radius R = AR, asymptotic temperature 
To = and fermion mass rrif = A The free energy is then rescaled as 

F = AF. Therefore, there exist a thermodynamic limit of N~‘^^^F, with N~‘^^^R, N^^^Tq, 
N^/^rrif approaching constant values when N —> 00 . In that limit the Thomas-Fermi equa¬ 
tion becomes exact. It is obvious that application of this model to stellar systems should 
work very well if the interactions among individual particles are negligible. 
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